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Abstract 

This  paper  presents  a  simplified  model  of  a  nickel  hydrogen  battery  based  on  a  first  approximation.  The  battery  is  assumed  uniform 
throughout.  The  reversible  potential  is  considered  primarily  due  to  one-electron  transfer  redox  reaction  of  nickel  hydroxide  and  nickel 
oxy hydroxide.  The  non-ideality  due  to  phase  reactions  is  characterized  by  the  two-parameter  activity  coefficients.  The  overcharge  process  is 
characterized  by  the  oxygen  reaction.  The  overpotentials  are  lumped  to  a  tunable  resistive  drop  to  fit  particular  battery  designs.  The  model 
is  implemented  in  the  Virtual  Test  Bed  environment,  and  the  characteristics  of  the  battery  are  simulated  and  in  good  agreement  with  the 
experimental  data  within  the  normal  operating  regime.  The  model  can  be  used  for  battery  dynamic  simulation  and  design  in  a  satellite  power 
system,  an  example  of  which  is  given. 

©  2004  Elsevier  B  .V.  All  rights  reserved. 
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1.  Introduction 

Complicated  electrochemical  processes  in  a  nickel  hydro¬ 
gen  battery  [1,2]  involve  more  than  eight  reactions  including 
main  reactions,  phase  reactions,  and  side  reactions.  While  a 
detailed  model  has  been  developed  [3,4],  it  is  currently  coded 
in  discipline- specific  languages  and  for  standalone  simula¬ 
tion  only.  Thus,  it  is  not  in  general  available  for  use  in  system 
level  simulation  and  design.  This  paper  presents  a  simpli¬ 
fied,  yet  still  physics-based,  model  that  can  be  rapidly  imple¬ 
mented  and  conveniently  used  for  the  battery  dynamic  sim¬ 
ulation,  design,  and  optimization  in  complicated  spacecraft 
power  systems. 

The  model  is  based  on  following  considerations.  For  typ¬ 
ical  applications  (e.g.,  in  spacecraft  power  systems  [5-7]), 
the  charge  and  discharge  of  the  Ni-FL  battery  are  in  such  a 
controlled  way  that  the  rate  of  charge  and  discharge  is  less 
than  0.5  C  to  ensure  a  reliable,  safe  and  long  lifetime  op¬ 
eration  [8,9].  Under  such  a  circumstance,  the  charge,  mass 
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and  heat  transport,  and  the  phase  equilibrium  can  be  as¬ 
sumed  to  complete  instantly.  Thus,  the  kinetics  of  reactions 
is  less  likely  a  limiting  factor  to  the  external  characteristics. 
This  leads  to  the  assumption  that  the  boundaries  between 
the  electrodes  and  the  separator,  and  those  between  the  gas 
phase  and  solid  phase  in  the  battery  vanish,  and  the  electro¬ 
chemical  and  electrothermal  processes  take  place  uniformly 
throughout  the  battery.  This  assumption  greatly  simplifies  the 
modeling  process  while  observing  the  major  features  of  the 
battery. 

The  present  model  is  implemented  in  the  Virtual  Test  Bed 
software  environment  [10],  in  which  the  battery  is  consid¬ 
ered  as  an  interdisciplinary  (chemical,  electrical  and  thermal) 
system  governed  by  physics  laws  under  certain  initial  and 
boundary  conditions.  Given  a  time,  the  battery  interacts  with 
surroundings  in  two  ways:  it  delivers  electrical  power  to  the 
load  at  its  electrical  terminals  during  discharge  by  convert¬ 
ing  chemical  energy  into  electric  one  (or  the  reverse  during 
charge  process),  and  it  absorbs  or  releases  heat  to  ambient 
at  its  external  surface.  Thus,  the  battery  voltage  is  a  result  of 
coupled  electrochemical,  chemical-thermal  and  electrother¬ 
mal  processes. 
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Nomenclature 

A 

electrode  area  (m2) 

a 

specific  area  (m2  m-3) 

bj(t  - 

h)  time  history  for  the  through  variable  at  termi- 
nal  j 

cb 

battery  capacity  (A  h) 

CP 

average  specific  heat  (Jkg-1  K-1) 

ENt 

total  reversible  potential  (V) 

Eo2 

oxygen  reaction  potential  (V) 

h 

simulation  time  step  (s) 

F 

Faraday  constant  (96,485  C  mol-1) 

gj,k 

admittance  of  a  through  variable  at  terminal  j 
with  respect  to  an  across  variable  at  terminal  k 

L+ 

electrode  thickness  (m) 

m 

battery  cell  mass  (kg) 

iNi 

nickel  reaction  current  (A) 

io2 

oxygen  reaction  current  (A) 

io2,o 

oxygen  exchange  current  density  (Am-2) 

h 

thermal  current  (W  K-1) 

«h2 

amount  of  hydrogen  (mol) 

p 

across  variable  for  the  pressure  terminal  (Pa) 

p 

pressure  (Pa) 

(2  max 

maximum  charge  stored  in  the  battery  (C) 

R 

gas  constant  (8.3 14  J  mol-1  K-1) 

Ri 

internal  resistance  (£2) 

T 

battery  bulk  temperature,  also  across  variable 
for  thermal  terminal  (K) 

t 

independent  variable  time  (s) 

Vg 

volume  of  battery  (m3) 

V 

battery  terminal  voltage  (V) 

Ca,  Vc 

across  variables  for  anode  and  cathode  (V) 

X 

state  of  discharge 

ao2 

oxygen  reaction  transfer  number  (0.75) 

The  solution  to  the  coupled  processes  is  obtained  by  means 
of  the  resistive-companion  (RC)  method  [11],  which  requires 
the  governing  equations  to  be  time-discretized  to  a  standard 
matrix  formulation  in  terms  of  across  and  through  variables. 
The  resulted  matrix  equations  are  then  simultaneously  solved 
through  the  VTB  time-domain  simulation  using  the  nodal 
analysis  approach. 


2.  Model  description 

2.7.  Electrochemical  process 

The  deviation  of  the  reversible  potential  of  the  nickel  elec¬ 
trode  from  the  ideal  Nernst  equation,  or  its  non-ideality,  was 
interpreted  [12]  as  due  to  intercalation  of  multi-phases,  in¬ 
cluding  a-  and  (3-phase  of  nickel  hydroxide,  and  (3-  and  y- 
phase  of  nickel  oxy hydroxide.  When  discharging,  it  yields 


higher  capacity  [4]  than  that  obtained  from  a  single-phase 
nickel  oxyhydroxide.  As  a  first  approximation,  the  overall  po¬ 
tential  can  be  characterized  as  if  a  usual  one-electron  transfer 
redox  reaction  occurs  between  nickel  hydroxide  and  oxy¬ 
hydroxide,  but  with  the  incorporation  of  the  effects  of  both 
non-ideality  due  to  multi-species  of  active  nickel  materials 
and  the  potential  due  to  hydrogen  reaction  at  the  platinum 
electrode. 

To  characterize  the  overcharge  process,  the  oxygen  reduc¬ 
tion  is  considered.  It  was  believed  [4]  that  the  water  oxida¬ 
tion/reduction  (or  oxygen  evolution)  at  the  nickel  electrode 
and  the  oxygen  reduction  at  the  platinum  electrode  plays  a 
dominant  role  in  an  overcharge  process,  causing  a  large  por¬ 
tion  of  the  current  wasted.  For  normal  charge  and  discharge, 
the  oxygen  reduction  yields  a  negligible  current. 

Based  on  the  above  descriptions,  the  electrochemical  reac¬ 
tions  that  characterize  the  energy  conversion  between  chemi¬ 
cal  and  electrical  forms  in  the  Ni-H2  battery  are  given  below. 

At  the  nickel  electrode  are  the  nickel  reduction/oxidation 
and  the  water  reduction/oxidation: 

charge 

Ni00H  +  H20  +  e-  <=►  Ni(OH)2  +  OH-,  (1) 

discharge 

1 

-02  +  H20  +  2e~  <=►  20H- .  (2) 

2 

At  the  platinum  electrode  are  the  hydrogen  reaction  and 
the  oxygen  reduction: 

\  charge 

-H2  +  OH"  <=>  H20  +  e-,  (3) 

^  discharge 

±02  +  H20  +  2e-  -*  20H”.  (4) 

The  overall  reversible  potential  can  be  described  as: 

En,(T.  x)  =  E°n.(T )  +  ^  In 

+  E-[Ao(T)C(x)  +  Bo(T)D(x)l  (5) 

2F 


1  - 


The  third  and  the  last  term  in  (5)  are  the  characteristics  of 
the  non-ideality,  and, 


A0(T)  =  a0  +  a\T,  (6) 

B0(T)  =  b0  +  biZ  (7) 

C(x )  =  co  +  c\x,  (8) 

D(x)  =  do  +  d\x  +  d2X2,  (9) 


eod  +  e\ T, 
6o  c  +  ci  T, 


for  discharge 
for  charge 


(10) 


where  <2;,  bu  c/,  and  d[  are  constants  for  phase  activities  of 
nickel  species,  which  may  vary  from  one  battery  to  another 
due  to  different  contents  and  pore  structures.  eoc,  co d  and  ci 
are  also  constants.  The  difference  between  eoc  and  eod  in  (10) 
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is  used  to  characterize  the  hysteresis  of  proton  intercalation 
in  nickel  hydroxide  electrode  [13,14]. 

Under  the  assumption  of  a  uniform  and  non-kinetic  limited 
reaction,  the  current  from  oxygen  reaction  can  be  calculated 
according  to  the  Butler- Volmer  equation, 


io2(v,  T )  =  AL+aio,o2(T ) 


x  exp 


4(1  -  ao2)—(v  -  Eo2(T )) 
K1 


(ii) 


where  the  exchange  current  density  is  of  the  Arrhenius  type 
[3,4],  and, 


E0i(T)  =  1.730-0.001687; 


(12) 


aio2,o(T)  =  10.0  exp 


-12,025 


(13) 


Eq.  (12)  is  the  simplified  equilibrium  potential  of  the  oxy¬ 
gen  reaction  depending  upon  the  temperature  only.  The  nu¬ 
merical  parameters  appeared  in  (12)  and  (13)  are  calculated 
based  on  the  data  in  references  [3,4]  that  are  specifically  ap¬ 
plied  to  the  oxygen  reaction  in  Ni-hU  batteries. 

The  state  of  discharge  v  can  be  conveniently  related  to  the 
nickel  reaction  current  as: 


dx  i  Nj 

dt  (2  max 


(14) 


The  battery  terminal  voltage  can  be  related  to  the  current: 


v  =  E^  +  Rtf  Nr  (15) 

The  fitting  parameter  Ri  accounts  for  the  overall  internal 
loss  due  to  ohmic,  surface  kinetic  and  concentration  limita¬ 
tions.  Since  the  internal  loss  in  the  Ni-H2  battery  is  gener¬ 
ally  small  compared  to  most  of  other  major  secondary  bat¬ 
teries,  we  simply  use  a  small  constant  to  approximate.  Eqs. 
(4)— (15)  complete  the  descriptions  of  the  electrochemical 
process. 


2.2.  Thermal  process 


To  simplify,  we  consider  the  heat  generations  due  to  the 
enthalpy  of  reactions  [15]  and  the  heat  transported  at  the  ther¬ 
mal  terminal  only,  whereas  the  heat  due  to  isometric  process, 
enthalpy  of  mixing,  phase  change  and  change  of  specific  heat 
are  ignored  without  introducing  significant  error.  The  rate  of 
temperature  change  of  the  battery  is  governed  by  the  energy 
balance  equation: 


d  T 


9ENi  \ 


cDm —  =  iN  \  v  —  Em.  +  T 
p  dt  '  V  '  dT  ) 


+  io2  (  v  ~  Eo2  +  T  )  +  ij  T. 


(16) 


The  first  two  terms  on  the  right-hand- side  are  owing  to 
the  enthalpy  of  reactions,  which  can  be  identified  as  (1) 
polarization  heat  including  /^(n  —  E^)  and  io2(v  —  Eq2), 


which  are  irreversible  electrothermal  processes;  and  (2)  en- 
tropic  heat  including  i^iT{dE^i/dT)  and  io2T(dEo2/dT ), 
which  are  reversible  chemical-thermal  processes.  The  last 
term  on  the  right-hand  side  is  the  heat  transported  to  (or 
from)  the  surroundings  through  the  battery  surface  or  the 
thermal  terminal,  which  is  expressed  as  the  product  of  the 
through  variable  ij  (thermal  current)  and  the  across  variable  T 
(temperature). 


2.3.  Pressure  modeling 


The  hydrogen  pressure  in  the  battery  is  an  indicator  for 
the  state  of  discharge  (hydrogen  content)  under  the  condition 
of  regulated  charge/discharge.  Thus,  a  built-in  sensor  is  used 
to  detect  the  pressure  [5],  so  that  a  VIP  (voltage/pressure) 
control  can  be  implemented  to  prevent  the  battery  from 
overcharge.  Assuming  the  ideal  gas  law  applies,  we  then 
have: 


(17) 


d^H2  _  iNi 
dt  ~  2 F’ 

3600Cb 

«h2I?=o  =  2/7  (1 


x\t=o) 


(18) 

(19) 


The  partial  pressure  of  oxygen  is  ignored  in  this  calcula¬ 
tion.  This  is  based  upon  the  assumption  that  the  reactions  in 
the  battery  are  uniform  throughout,  and  the  charge/discharge 
is  not  so  fast  that  any  appreciable  oxygen  produced  will  be 
immediately  removed  via  reaction  (4). 


3.  VTB  implementation 

The  VTB  implementation  of  the  model  described  by 
Eqs.  (5)— (19)  involves  (1)  defining  appropriate  terminals 
for  power  or  signal  coupling;  (2)  deriving  the  terminal 
equations  in  their  standard  RC  format  for  each  terminal; 
and  (3)  developing  C++  code  for  computer  execution.  The 
time-domain  solution  for  coupled  physics  processes  of  the 
battery  can  then  be  obtained  through  VTB  simulation  for 
any  given  circuit  configurations  to  which  the  battery  is 
connected. 

The  physical  electric  terminals  of  the  Ni-EU  battery  are 
its  anode  and  cathode  that  can  be  connected  to  an  external 
circuit  for  electrical  power  coupling.  Naturally,  the  current 
for  the  through  variable  and  the  potential  for  the  across  vari¬ 
able  are  defined.  The  physical  thermal  terminal  is  the  exter¬ 
nal  surface  through  which  the  heat  can  transfer  to  or  from 
the  surroundings.  By  definition,  the  thermal  terminal,  or  the 
battery  surface  can  transport  heat  power  by  means  of  con¬ 
duction,  or  radiation,  or  convection.  The  pressure  terminal 
is  a  signal  terminal.  The  signal  from  a  built-in  pressure  sen¬ 
sor  is  coupled  to  an  external  circuit  through  this  signal  port 
for  monitoring  or  control  purposes.  Thus,  there  is  no  actual 
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Fig.  1.  VTB  schematic  view  of  the  N1-H2  battery. 


mass  or  energy  flow  at  this  terminal.  This  is  the  major  dif¬ 
ference  between  a  signal  terminal  and  a  natural  terminal, 
where  the  power  is  transported,  such  as  the  anode,  cathode, 
and  thermal  terminal  defined  here.  VTB  implementation  of 
the  pressure  signal  terminal  is  to  convert  the  pressure  into 
a  voltage  signal  while  allowing  an  infinitesimal  current  at 
the  terminal.  This  is  shown  in  Appendix  A.4.  The  defini¬ 
tion  of  the  terminals  and  terminal  variables,  together  with 
the  battery  icon  in  the  VTB  schematic  view,  are  shown  in 
Fig.  1. 

Observing  that  the  electrochemical  and  electrothermal 
processes  are  coupled,  the  through  variables  are  generally 
functions  of  all  possible  across  variables.  Specifically,  we 
expect  that, 


( 4(0  ^ 

^  £a,a 

8  a,c 

ga,T 

0  ^ 

/  Pa(0\ 

4(0 

£c,a 

Sc.c 

<?c,T 

0 

Pc(0 

4(0 

£T,a 

<?T,c 

gT.T 

0 

Tit) 

VgP.a 

gP.c 

gP,T 

gp,p/ 

V  P(t) ) 

/  K(t  -h)\ 


bj(t  —  h) 

\bP(t-h)J 


The  four- terminal  RC  model  for  the  Ni-H2  battery  given 
by  (22)  is  coded  in  C++  language,  and  then  compiled  into  a 
model  object.  Within  the  VTB,  the  battery  model,  appeared 
as  that  in  Fig.  1,  can  be  easily  inserted  into  any  circuit  for 
simulation.  For  the  vector  Eq.  (22)  to  accurately  represent 
the  characteristics  of  the  battery,  the  calculation  of  admit¬ 
tance  matrix  G(t  —  h )  and  the  time  history  vector  B(t  —  h) 
is  critical.  During  simulation,  the  VTB  solver  seeks  a  solu¬ 
tion  to  the  nodal  variables  for  the  entire  system  for  each  time 
step.  The  results  as  functions  of  time  can  be  displayed  on 
the  oscilloscope-like  screens,  or  can  be  saved  as  data  files  for 
later  use.  Eq.  (22)  can  yield  satisfactory  results  provided  that 
the  simulation  time  step  is  sufficiently  small. 


4.  Battery  characteristics 

4.1.  Parameter  setup 


ij(t )  =  -  bj(t  -  h ),  j,  k  =  va,  vc,  T,  P, 

k 

(20) 


where, 


,  and  bj(t  —  h ) 

t—h 

=  -ijit  -  h)  +  j,kk(t  -h). 

k 

In  above  equations,  j  signifies  for  subscripts  only,  whereas 
k  is  used  to  denote  both  subscripts  and  across  variables  of 
na,  vc,  T  and  P.  In  matrix  notation,  Eq.  (20)  is  of  the  following 
form: 

I(t)  =  G(t  -  h)V(t)  +  B(t  -  h ),  (21) 

where  /(0  =  {4(0Uxi>  V(t)  =  {k(t)} 4xl,  G(t  -  h)  = 
{gj,k} 4x4’  and  B(t  -  h)  =  {bj(t  -  h)}4xl,  j,  k  — 
t>a,  Vc,  T,  P. 

The  result  for  the  RC  equations  of  the  electrical  terminals 
(anode  and  cathode),  thermal  terminal  and  pressure  signal 
terminal,  based  on  the  mathematical  Eqs.  (5)— (19),  are  sum¬ 
marized  in  Eq.  (22)  (the  detailed  derivation  is  provided  in 
Appendix). 


Sj,k  — 


<2i 

dk 


We  first  demonstrate  the  Ni-H2  battery  static  characteris¬ 
tics  using  the  circuit  shown  in  Fig.  2,  where  the  battery  model 
is  subject  to  a  charge-overcharge-self-discharge-discharge 
test.  The  electrodes  are  connected  to  a  current  load  that  is  ca¬ 
pable  of  sourcing  or  sinking  current.  A  heat  sink  with  a  large 
heat  transfer  coefficient  is  connected  to  the  thermal  terminal 
to  maintain  the  battery  temperature  in  nearly  a  constant.  The 
heat  sink  may  transport  heat  from  the  battery  to  the  ambi¬ 
ent,  or  vice  versa.  The  pressure  signal  terminal  is  left  open 
at  the  moment.  The  model  represents  a  single  battery  cell 
consisting  of  a  stack  of  24  modules  connected  in  parallel 

NiH2  Battery 


Heat  sink 


Ambient  Temp. 


Fig.  2.  VTB  simulation  circuit  for  the  M-H2  battery  characteristics. 
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Table  1 


Parameter  setup  for  a  single  Ni-FU  battery  cell  [4] 


Parameter 

Value 

Unit 

Cell 

Capacity 

30 

Ah 

Volume 

4.65  x  1CT4 

m3 

Mass 

1.2 

kg 

Resistance 

0.01 

£2 

Surface  area 

0.06 

2 

m 

Specific  heat 

1000.0 

«— 1 

CfQ 

1 

R 

1 

No.  of  modules  in  series  24 

Module 

Electrode  area 

0.005 

2 

m 

Electrode  thickness 

3.6  X  1(T4 

m 

Table  2 

Parameters  for  multi-phase  activity 

A0(T) 

ao  =  1 1 .50,  a\  =  —0.023 1 

Bo(T) 

b0  =  32.00,  bi=  0.060 

C(x) 

cq  =  4.590,  ci  =  15.80 

D(x) 

d0  =  1.800,  d\  =  —5.680,  J2  =0.350 

E%(T ) 

cqc  =  1 .609,  cq d  -  1  -470,  ci  =  0.00062 

and  contained  in  a  cylindrical  pressure  vessel.  Each  module, 
characterized  by  Eq.  (22),  has  a  nickel  electrode,  platinum 
electrode,  Zircar  separator  and  polypropylene  gas  screen  (for 
detailed  constructions,  please  see  references  [2, 4, 8, 9]).  Other 
physical  parameters  are  listed  in  Table  1 .  The  parameters  for 
multi-phase  activity  are  given  in  Table  2.  In  the  simulation, 
the  initial  temperature  is  set  to  the  ambient  temperature,  and 
the  initial  state  of  discharge  is  set  to  100%. 

4.2.  Temperature -dependent  characteristics 

The  simulation  was  first  conducted  for  a  cycle  of  16  h 
charge  at  a  C/10  rate  (3  A),  lh  self-discharge,  and  then  a 
discharge  at  a  C/2  rate  (15  A)  till  it  is  depleted.  The  total 
runtime  is  19  h.  Shown  in  Fig.  3  is  the  load  current  profile. 


4 

1  1  1  1 

2 

- 

0 

- 

-2 

m 

=  -4 

~ 

fc 

5  -6 

■ 

b 

2  -8 

■ 

- 

CQ 

CQ 

-10 

■ 

-12 

- 

-14 

■ 

- 

-16 

,i - . - 1 - 1 - 

0  5  10  15  20 


Time  (hour) 

Fig.  3.  The  load  current  profile. 


Fig.  4.  Cell  voltages  for  a  cycle. 


The  cycle  was  repeated  for  three  ambient  temperatures  —273, 
283  and  293  K. 

The  battery  cell  voltage  profiles  (lines)  are  shown  in  Fig.  4 
and  compared  against  the  TRW  30  A  h  cell  data  (points)  [4] . 
As  can  be  seen,  the  model  fits  the  experimental  data  quite 
well  for  the  normal  charging  (first  10  h)  at  all  three  tempera¬ 
tures.  The  model  is  also  in  good  agreement  with  the  data  for 
self-discharge  and  discharge  above  1.15V(16-19h).  During 
overcharges,  the  model  yields  constant  voltages,  while  the 
experimental  data  show  slow  ramps.  This  is  primarily  due  to 
the  temperature-dependent  activities  of  the  multi- species  of 
nickel  materials  and  possibly  the  non-linear  overpotentials, 
which  are  not  well  characterized  by  the  model.  However,  the 
ultimate  overcharge  voltages  from  the  model  are  matched 
to  those  from  the  experiments.  In  addition,  the  simulated 
voltage  behaviors  are  consistent  with  the  data  over  the  en¬ 
tire  cycle.  That  is,  lower  temperatures  yield  higher  voltages 
and  also  longer  discharge  life,  suggesting  that  the  battery  is 
preferably  operated  at  low  temperatures.  The  discrepancy  be¬ 
tween  the  model  and  the  data  is  also  seen  when  discharging 
below  1.15  V,  where  the  experimental  data  show  larger  dis¬ 
charge  capacities.  This  is  due  to  the  appearance  of  y-phase 
of  nickel  oxyhydroxide  [4],  which  will  be  discussed  further 
later.  Nonetheless,  it  can  be  concluded  that  the  model  matches 
the  experimental  data  very  well  within  the  normal  operation 
range. 

Figs.  5  and  6  show  the  currents  from  the  nickel  reaction 
and  the  oxygen  reaction  respectively  for  the  same  cycles.  At 
any  moment,  the  sum  of  the  nickel  current  and  the  oxygen 
current  is  equal  to  the  total  current.  These  figures  clearly 
show  how  a  dominant  nickel  reaction  in  normal  charge  pro¬ 
cess  transits  to  a  dominant  oxygen  reaction  in  an  overcharge 
process.  During  overcharges,  the  majority  of  a  charge  current 
is  wasted  for  the  oxygen  reaction  since  the  energy  is  not  used 
for  oxidizing  nickel  hydroxide.  During  a  self-discharge  pro¬ 
cess,  the  source  of  the  oxygen  current  is  the  reduction  of  the 
nickel  oxyhydroxide,  thus  yielding  nickel  and  oxygen  cur- 
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Fig.  5.  The  nickel  reaction  currents  for  a  cycle. 


rents  of  equal  amount  but  opposite  directions.  For  a  normal 
discharge,  the  oxygen  current  is  usually  minimal  although  it 
can  be  significant  at  high  temperatures  (see  curves  for  293  K). 

The  state  of  discharge,  shown  in  Fig.  7,  varies  almost  lin¬ 
early  for  a  constant  charge  or  discharge  rate.  During  an  over¬ 
charge,  the  state  of  discharge  maintains  nearly  a  constant 
value.  The  value  is  higher  for  a  higher  temperature  because 
more  current  goes  to  the  oxygen  reaction.  Consequently,  a 
lower  overvoltage  plateau  results  (see  Fig.  4).  The  pressure 
profiles  shown  in  Fig.  8  are  complementarily  similar  to  the 
states  of  discharge.  This  is  in  fact  due  to  the  stored  hydrogen 
representing  the  charge  state.  However,  since  the  pressure  is 
not  only  proportional  to  the  amount  of  hydrogen,  but  also  to 
the  temperature,  it  does  not  accurately  represent  the  state  of 
discharge.  For  example,  the  pressures  are  higher  for  higher 
temperatures  for  almost  the  entire  cycle  rather  than  lower  to 
be  consistent  with  the  charge  state  behaviors,  except  at  the 
joining  area  of  normal  charge  and  overcharge  where  a  signif¬ 
icant  current  loss  to  the  oxygen  reaction  causes  the  pressure 
for  293  K  to  be  lower  than  the  pressure  for  283  K.  However, 


Fig.  7.  The  state  of  discharge  for  a  cycle. 


this  does  not  generally  affect  the  use  of  the  pressure  detec¬ 
tion  for  overcharge  protection.  Thus,  our  simplified  pressure 
modeling  can  still  serve  the  purpose. 

The  temperatures  of  the  battery  for  the  cycle  test  are 
recorded  in  Fig.  9.  Notice  that  the  TRW  30  Ah  cell  has  an 
effective  heat  transfer  area  of  0.06  m2.  For  maintaining  a 
constant  battery  temperature  during  the  test,  the  cooling 
coefficient  can  be  as  large  as  60  Wm~2  K_1  [4],  which  is 
reasonable  for  a  forced  convection  cooling.  In  our  simulation, 
the  rate  of  heat  transfer  was  set  equivalently  to  3.6  WK_1, 
so  the  conditions  for  cooling  are  comparable  to  the  experi¬ 
mental  condition.  Fig.  9  shows  the  battery  cell  temperatures 
are  slightly  higher  than  the  ambient  during  charging  for 
each  test,  but  with  small  variations.  During  discharge,  the 
temperature  increases  significantly  due  to  two  reasons: 
(1)  the  overpotential  heat  due  to  a  large  discharge  current, 
and  (2)  the  reversible  heat  release  during  electrochemical 
reactions.  Clearly  the  heat  generation  rate  exceeds  the 
surface  transfer  rate  and  therefore  causing  the  temperature 


Fig.  6.  The  oxygen  reaction  currents  for  a  cycle. 


Fig.  8.  The  battery  pressures  for  a  cycle. 
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Fig.  9.  The  battery  cell  temperature  variations  for  a  cycle. 


to  rise.  Because  of  a  large  cooling  coefficient  and  a  small 
cell  mass,  a  thermal  equilibrium  between  the  battery  and  the 
ambient  was  quickly  reached.  This  is  shown  by  the  fact  that 
the  cell  temperature  reaches  a  constant  value  soon  after  the 
battery  transits  from  one  steady- state  mode  to  another. 

4.3.  Rate-dependent  characteristics 

Figs.  10-13  show  the  battery  characteristics  under  three 
charge/discharge  rates:  0.5,  1.0  and  2.0  C.  The  ambient  tem¬ 
perature  is  283  K.  For  the  three  rates,  the  battery  was  charged 
for  2.2, 1 . 1  and  0.55  h  respectively,  and  then  discharged  for  2, 
1  and  0.5  h,  respectively.  In  Fig.  10,  both  the  voltages  (thick 
lines)  and  the  reversible  potentials  (thin  lines)  are  shown. 
For  discharging  at  high  rates,  the  battery  voltages  are  consid¬ 
erably  lower  than  its  reversible  potentials  due  to  significant 
internal  loss  that  is  typical  to  all  batteries.  For  overcharge, 
the  holding  potentials  are  higher  for  higher  rates  due  to  a 
lower  temperature  (reversible  heat  effect).  Figs.  11  and  12 


Fig.  10.  The  cell  voltages  (thick  lines)  and  the  reversible  potentials  (thin 
lines). 


Fig.  11.  The  state  of  discharge  for  different  rates. 


Fig.  12.  The  cell  pressures  for  different  rates. 


Fig.  13.  The  battery  temperatures  for  different  rates. 
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are  the  correspondent  states  of  discharge  and  pressures  for 
the  three  rates,  showing  similar  characteristics  to  those  in 
previous  simulations. 

Fig.  13  shows  the  cell  temperature  for  the  three  rates.  It 
is  lower  than  the  ambient  during  charging  processes,  higher 
than  the  ambient  in  discharging  processes.  Apparently,  the 
reversal  of  the  entropic  heat  in  charging  process  causes  the 
temperature  to  drop.  During  discharging,  the  release  of  en¬ 
tropic  heat  and  irreversible  heat  result  in  higher  temperatures. 
It  is  worth  to  remark  that  the  temperature  drop  during  charg¬ 
ing  was  also  noticed  in  reference  [4].  Our  next  work  is  to 
improve  the  thermal  modeling  once  the  data  are  available  for 
characterization. 


5.  Battery  dynamics  in  a  system 

5.7.  System  description 

The  dynamic  behavior  of  the  Ni-H2  battery  in  a  simpli¬ 
fied  spacecraft  electric  power  system,  established  in  VTB  as 
shown  in  Fig.  14,  is  simulated.  The  spacecraft  is  a  low-earth- 
orbit  (LEO)  satellite  having  an  average  altitude  of  798  km, 
an  inclination  angle  of  60°,  and  an  orbit  of  about  100  min,  of 
which  35  min  are  the  shadow  time. 

The  energy  source  for  the  system  is  a  solar  array  (SA). 
The  Ni-H2  battery  array  stores  energy  at  the  sun-view  time, 
and  provides  power  to  the  load  at  the  earth-shadow  time. 
In  order  for  the  battery  to  operate  properly,  a  constant- 
current/constant-voltage/pressure-limited  algorithm  is  em¬ 
ployed  for  charge  control  and  protection.  The  algorithm  is 
implemented  using  the  following  components:  a  step-down 
(buck)  converter,  a  PID  controller,  sensors,  and  diodes.  The 
details  of  these  component  models  can  be  found  in  VTB  web¬ 
site  [16],  will  not  be  described  here.  The  signals  of  the  battery 
voltage  and  current  are  sent  to  the  controller  for  appropriate 
execution  of  the  algorithm.  In  addition,  the  pressure  signal 
feeds  the  controller  for  overcharge  protection.  Heat  sinks  to 
both  the  solar  array  and  the  battery  transport  the  heat  to  am¬ 
bient.  For  simplicity,  a  single  load  -  transmitter  -  is  used  in 
the  system. 


Table  3 


Parameter  setup  for  the  spacecraft  power  system 


Parameter 

Specification 

Spacecraft  orbit 

Altitude  798  km,  inclination  60°,  orbit  period 

100  min 

Solar  array 

100S  x  3 OP  cells.  Cell  Voc  =  0.615  V,  area 

0.01  m2,  series  resistance  10m£2 

Battery  array 

24S  x  8P  cells.  The  cell  parameters  are  given  in 
Table  1 

Controller 

Constant-current  32  A,  constant- voltage  36  V, 
pressure  limit  2.6  x  106  Pa 

Transmitter 

3.25  kW  active,  0.25  kW  sleeping.  Duty  ratio 

50%,  period  17.5  min 

The  system  configuration  is  not  necessarily  optimized. 
However,  the  solar  array  size  is  chosen  such  that,  without 
shunt  regulation,  it  delivers  a  power  to  the  load  and  the  bat¬ 
tery  at  a  high  efficiency  (1 1.5%  at  a  high  power  phase).  The 
battery  is  sized  by  considering  (1)  it  stores  and  delivers  a 
nearly  equal  amount  of  energy  in  each  orbit  without  experi¬ 
encing  a  deep  discharge  hazard;  (2)  the  capacity  of  the  battery 
is  designed  to  meet  the  end  of  life  (EOL)  requirement  based 
on  the  Space  Systems  Loral  performance-based  battery  de¬ 
sign  life  model  [8,9].  That  is,  with  a  temperature  of  283  K 
and  a  maximum  depth  of  discharge  of  35%,  the  model  yields 
a  battery  life  of  about  6.5  years.  The  parameter  setup  for  the 
major  components  of  the  system  is  listed  in  Table  3.  These 
numbers  are  chosen  to  meet  the  above  criteria  for  demonstra¬ 
tion  and  observation  of  the  battery  dynamic  behaviors  in  a 
simplified  system  environment. 

5.2.  Simulation  results 

The  total  simulation  time  is  233.3  min,  about  2.33  or¬ 
bits.  The  results  are  shown  by  Figs.  15-23.  Fig.  15  shows 
the  orbit  solar  insolation  received  by  the  solar  array  while 
it  is  orbiting.  The  insolation  level  is  1357Wm-2  when  the 
array  sees  the  sun,  and  it  drops  to  zero  at  the  shadow  pe¬ 
riod.  The  load  power  drawn  by  the  transmitter  is  shown 
in  Fig.  16.  Notice  that  the  transmitter  was  operated  at  the 
specified  power  levels  and  cycle  rate  regardless  of  the  in¬ 
solation  levels  at  the  solar  array.  Obviously,  at  the  shadow 


Fig.  14.  Spacecraft  power  system  using  N1-H2  battery  as  an  energy  storage  device. 
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Fig.  15.  Insolation  received  by  the  solar  array  as  a  function  of  orbiting  time. 


Fig.  17.  The  solar  array  output  power  for  2.33  orbits. 


time,  the  transmitter  power  was  provided  by  the  battery 
array. 

Fig.  17  shows  the  solar  array  power.  At  the  sun- view  time, 
the  solar  array  converts  the  solar  energy  and  sends  power 
to  the  load  as  well  as  to  the  battery.  Due  to  a  pulsed  load, 
the  output  power  is  also  pulsating.  For  a  high-power  load  of 
3.25  kW,  the  solar  array  power  is  4.64  kW  with  an  efficiency 
of  11.5%  and  a  voltage  of  59.5  V.  For  a  low-power  load  of 
0.25  kW,  the  array  power  is  about  1 .42  kW  with  an  efficiency 
of  3.5%  and  a  voltage  of  55.4  V.  At  the  earth-shadow  time, 
the  solar  array  power  is  zero,  leaving  the  battery  to  supply 
the  power  to  the  load.  The  details  of  the  solar  array  character¬ 
istics,  including  voltage,  current,  efficiency  and  temperature, 
are  available  from  the  simulation,  but  will  not  be  covered  here 
since  it  is  not  a  focus  of  this  paper. 

The  battery  behaviors  in  orbit  are  shown  in  Figs.  18-23. 
At  the  sun- view  time,  the  battery  is  charging,  having  an  in¬ 
creasing  voltage  from  34.3  to  34.7  V  (Fig.  18)  and  a  constant 
current  of  32  A  (Fig.  19),  a  rate  of  C/7.5.  The  positive  power 
(about  1 . 1  kW,  Fig.  20)  indicates  that  the  battery  is  absorb- 


Fig.  18.  The  battery  array  voltage  for  about  2.33  orbits. 
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Fig.  16.  The  transmitter  power  profile  as  a  function  of  time. 


Fig.  19.  The  battery  array  current  for  about  2.33  orbits. 
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Fig.  20.  The  battery  array  power  for  about  2.33  orbits. 


Fig.  22.  The  battery  pressure  for  about  2.33  orbits. 


in g  power.  The  power  relations  among  the  source  power  in 
Fig.  17,  the  load  power  in  Fig.  16,  and  the  stored  power  in 
Fig.  20  reveal  that  there  is  a  power  loss  about  0.29  kW  at  a 
high-power  load,  and  a  loss  of  0.07  kW  at  a  low-power  load, 
which  are  caused  by  the  resistances  of  circuit  components. 

At  the  earth-shadow  time,  the  battery  takes  the  duty  to  de¬ 
liver  power  to  the  load  at  a  voltage  around  31V.  The  current 
has  large  ripples  due  to  a  pulsed  power  load.  For  a  high-power 
load  of  3.25  kW,  the  battery  discharge  current  is  130  A,  equiv¬ 
alent  to  a  rate  of  approximately  C/2;  while  for  a  low-power 
load  of  0.25  kW,  the  averaged  discharge  current  is  12.8  A, 
a  rate  of  about  C/19.  The  battery  output  power  (negative) 
is  4  kW  for  a  high-power  phase  and  averaged  0.4  kW  for  a 
low-power  phase.  The  system  losses  are  0.75  kW  for  high- 
power  phase,  and  0.15  kW  for  a  low-power  phase.  These 
losses  are  primarily  due  to  large  battery  currents  while  oper¬ 
ating  at  a  low  bus  voltage,  as  the  system  is  not  an  optimized 
design. 

Figs.  21  and  22  show  the  state  of  discharge  and  the  pres¬ 
sure  respectively,  which  restore  their  initial  values  at  the 


end  of  each  orbit,  since  the  total  energy  stored  and  the  de¬ 
livered  in  one  orbit  are  nearly  balanced.  The  state  of  dis¬ 
charge  at  the  end  of  an  orbit  is  about  0.3,  equivalent  to 
a  depth  of  discharge  70%,  well  above  the  depth  of  dis¬ 
charge  35%  requirement  for  6.5  years  life  according  to  the 
Space  Systems  Loral  performance-based  battery  life  model 
[8,9]. 

The  battery  temperature  variation  during  the  course  of  or¬ 
biting  is  shown  in  Fig.  23.  The  thermal  environment  in  the 
LEO  is  a  result  of  solar  insolation,  earth  albedo,  earthshine, 
spacecraft  heat  generation,  and  space  heat  sink  [17].  In  this 
study,  an  equilibrium  state  for  the  thermal  environment  is 
assumed  for  simplicity,  which  possesses  a  temperature  of 
273  K.  For  the  Ni-H2  battery  to  have  a  better  cycle  life,  it  is 
important  to  control  the  battery  not  to  operate  in  high  temper¬ 
ature  regimes.  We  have  chosen  the  heat  dissipater  such  that 
the  temperature  of  the  battery  does  not  exceed  283  K  dur¬ 
ing  operation  to  satisfy  6.5  years  life  requirement.  As  shown 
in  Fig.  23,  the  battery  temperature  has  an  about  10  K  swing 
(270-280  K)  for  each  orbit.  The  temperature  decreases  for 


Fig.  21.  The  state  of  discharge  of  the  battery  for  about  2.33  orbits.  Fig.  23.  The  battery  temperature  for  about  2.33  orbits. 
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a  charging  mode.  This  is  again  due  to,  in  addition  to  the 
heat  transfer,  the  electrochemical  reactions  that  absorb  the 
reversible  heat.  During  a  discharge  mode,  heat  is  released 
from  the  both  reversible  and  irreversible  processes,  causing 
the  temperature  to  rise.  However,  for  a  low  discharge  current, 
the  rate  of  heat  dissipation  is  higher  than  that  of  generation, 
yielding  a  temporary  temperature  drop.  This  is  shown  by  the 
jagged  temperature  profile  when  the  spacecraft  is  in  the  earth 
shadow  time. 


6.  Discussion 

The  present  model  for  the  Ni-H2  battery  is  based  on  a 
first  approximation — uniform  reactions,  no  rate  and  concen¬ 
tration  limitations,  and  the  two-parameter  activity  coefficient 
characterization  of  phase  activities.  The  model  allows  to  pre¬ 
dict  the  battery  behaviors  and  to  perform  preliminary  system 
level  design  and  analysis  under  the  conditions  of  a  regulated 
charge/discharge  protocol.  However,  like  any  other  models, 
the  present  model  has  several  limitations. 

The  first  noticeable  limitation  is  that,  as  shown  in  Fig.  4, 
the  model  does  not  agree  with  the  data  for  the  discharge 
below  1.15  V.  The  discharge  plateau  appeared  near  1.15  V 
is  due  to  the  activity  of  the  y-phase  nickel  oxyhydroxide, 
for  which,  apparently,  the  two-parameter  activity  coefficient 
method  does  not  characterize  very  well.  The  charge  stored 
in  the  y-phase,  according  to  [4],  makes  up  about  30%  of 
the  total  active  material.  The  conversion  to  the  y-phase  oc¬ 
curs  mainly  in  an  overcharge  process  at  a  cost  of  decreas¬ 
ing  the  (3-phase  oxyhydroxide.  During  discharge,  two  major 
parallel  processes  compete:  one  is  the  conversion  from  the 
|3 -phase  oxyhydroxide  to  the  hydroxide,  the  other  is  from  the 
y-phase  oxyhydroxide  to  the  hydroxide,  and  both  are  strongly 
temperature-dependent  and  contribute  to  the  discharge  capac¬ 
ity.  On  the  other  hand,  charging  at  high  temperatures  does 
not  yield  a  complete  oxidation  of  the  hydroxide,  which  con¬ 
versely  affects  the  discharge  capacities.  This  is  shown  by  the 
data  in  Fig.  4  that  a  larger  discharge  capacity  results  at  a 
lower  temperature,  and  the  voltage  discrepancy  is  decreased 
at  high  temperatures.  Without  an  overcharge,  the  fraction  of 
the  y-phase  is  less  than  10%,  whereas  the  (3-phase  is  nearly 
90%.  Based  on  this  analysis,  we  can  conclude  that  the  present 
model  best  applies  to  the  case  in  which  the  battery  is  operated 
at  a  low  temperature  (between  273  and  283  K),  and  is  con¬ 
trolled  such  that  it  does  not  undergo  an  overcharge  process. 
Otherwise,  the  model  incurs  an  error  of  30%,  for  which  the 
lower  voltage  limit  (1.15  V)  corresponds  to  an  actual  depth 
of  discharge  of  about  30%.  It  is  possible  in  the  future  to  char¬ 
acterize  the  phase  activities  using  three-  or  four-parameter 
activity  coefficients  for  improving  the  model  behaviors  in 
overcharge  and  deep  discharge  regimes. 

The  thermal  modeling  mainly  considers  the  constant  over¬ 
potential  loss  and  the  reversible  heat,  which  best  applies  to 
the  case  of  a  controlled  charge/discharge  scenario  described 
in  Section  5.  The  bulk  temperature  is  usually  lower  than  that 


of  porous  electrodes  where  most  of  heat  is  generated.  Thus, 
the  model  is  prone  errors  if  overcharge,  deep  discharge,  and 
high-rate  discharge  occur.  The  error  may  also  result  from 
following  sources.  (1)  Several  thermal  effects  are  ignored  in 
the  model,  such  as  heat  due  to  isometric  process,  enthalpy  of 
mixing,  phase  change  and  change  of  specific  heat.  (2)  The 
non-linear  feature  of  overpotential  loss  is  not  included.  (3) 
The  two-parameter  activity  coefficient  method  may  also  re¬ 
sult  in  discrepancy  in  reversible  heat.  These  features  can  be 
included  and  improved  in  the  future  modeling  using  detailed 
thermodynamic  data,  which,  unfortunately,  are  not  available 
at  this  time. 

Finally,  the  model  does  not  include  capacity  fade  effect. 
However,  the  battery  array  can  be  sized  based  on  the  Space 
Systems  Loral  performance-based  battery  design  life  model 
[8,9]  using  the  present  model,  as  explained  in  Section  5.  To 
ensure  a  safe,  reliable  and  long  lifetime,  the  spacecraft  bat¬ 
tery  discharge  rate  is  controlled  not  to  exceed  C/2.0,  the  dis¬ 
charge  voltage  not  to  exceed  the  lower  limit  of  1.15  V,  and 
the  temperature  not  to  exceed  283  K.  The  system  design  can 
be  improved  if  a  practical  load  profile  is  given  by  utilizing  an 
appropriate  bus  voltage  level  and  a  regulation  method,  and 
by  optimizing  solar  array  and  battery  array  sizes.  Though  it 
is  not  optimized,  the  present  system  still  allows  performing 
preliminary  studies  of  the  dynamic  behaviors  of  the  Ni-H2 
battery  in  a  LEO  space  environment. 

7.  Conclusion 

The  Ni-H2  battery  model  based  on  the  thermodynamics 
of  one-electron  redox  reaction  between  nickel  hydroxide  and 
nickel  oxyhydroxide  is  established  in  the  VTB.  The  battery 
is  assumed  uniform  throughout,  and  it  is  not  subject  to  rate 
and  concentration  dependencies,  which  allows  great  simpli¬ 
fication  in  mathematical  treatment.  The  non-ideality  of  the 
reversible  potential  is  characterized  by  the  two-parameter  ac¬ 
tivity  coefficients.  The  overcharge  and  the  self-discharge  pro¬ 
cesses  are  characterized  by  the  oxygen  reaction.  The  details 
of  the  energy  conversion  processes  are  discussed  and  the  VTB 
implementation  of  the  model  is  presented.  The  model  allows 
for  the  battery  behaviors  to  be  studied  in  a  configurable  sys¬ 
tem  while  it  is  dynamically  interacting  with  the  system  and 
the  environment.  The  simulation  results  show  that  the  model 
in  general  has  a  good  agreement  with  the  experimental  re¬ 
sults  above  a  voltage  of  1 . 15  V,  below  which  the  discrepancy 
appears  if  the  battery  is  overcharged  prior  to  discharge.  The 
future  research  will  be  devoted  to  improving  the  character¬ 
ization  of  thermal  behaviors,  non-linear  overpotentials,  and 
the  phase  activities. 
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Appendix  A 

A.l.  Useful  expressions 


Since  the  analytical  expressions  for  the  potentials  of  the 
nickel  and  oxygen  reactions  are  given  (Eq.  (5)  and  (12)),  they, 
together  with  their  derivatives,  can  be  accurately  computed. 
Thus,  we  will  seek  expressions  for  the  admittance  matrix  el¬ 
ements  expressed  in  terms  of  potentials  and  their  derivatives, 
which  are  given  as: 

dpNj 

dx 

=  -f  ((1 1~  -  -  X-[Aq{T)C'(x)  +  B0(T)D'(x)U  , 

(Al) 


9£/v;  A  A',.  E 

- -  =  e\^ - - 

dT  T 

+  U?[A,0(T)C(x)  +  B'0(T)D(x)l 
d2  E  n-  R 

-jfr  =  j[A'0(T)C(x)  +  B'0(T)D(x)], 


(A2) 

(A3) 


3  xdT 


R 

F 


/ _ 

x(l 


—  -  l-[A0(T)C\x)  +  Bo(T)D'(x)] \ 

—  x)  2 

-Ua'0(T)C'(x)  +  B'0(T)D'(x)1  j 

(A4) 


d£’o, 

— =  -0.00168. 
dT 

A.  2.  Electrical  terminal 


(A5) 


The  terminal  variables  can  be  related  to  the  battery  current 
and  voltage  as 

4(0  =  -4(0  =  4v,(0  +  4>2(0,  (A6) 

Va(0  -  vc(t)  =  00  =  ENi{t)  +  RiiNiit)-  (A7) 

Since  the  state  of  discharge  is  not  a  terminal  variable,  it 
is  necessary  for  the  device  object  to  compute  v  during  each 
time  step.  To  do  so,  Eq.  (14)  is  discretized  in  the  time  domain 
using  Gear’s  2  [18], 

iNi(t)  =  ——^E[2>x(t)  -  4 x(t  -h)  +  x(t  -  2ft)].  (A8) 


Combining  Eq.  (A7)  and  (A8),  we  have. 

t>a(0  -  Veit )  -  £0,  (0 

Bint 

=  ——^[},x{t)  -  4 x(t  -h)  +  x(t  -  2 h)].  (A9) 

2  h 

Eq.  (A9)  can  be  solved  for  x(t)  numerically  for  given  ter¬ 
minal  across  variables. 

The  nickel  current  component  can  be  derived  from  Eqs. 
(5)  and  (A7).  This  is, 

Z,v,(4  =  T: gNj,kk(t)  -  bNi(t  -h),  k  =  Da,  vc,  T  P, 

(A10) 


where, 
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.  gNt, P  =  0 


Tint  9T 


(All) 


t—h 


bNi(t  -h)=  ~h)  +  f2  8Nhkk(t  -  h), 


k  =  Da,  vc ,  T,  P. 


(A  12) 


For  the  oxygen  current  component,  since  it  is  an  explicit 
function  of  the  terminal  variables,  it  is  straightforward  to  find 
the  standard  expression,  as  given  below. 

4)2 (0  =  T:  gQ2,kk(t )  -  bo2(t  -h),  k  =  va,  ve,  T  T, 

(A  13) 


where, 


^02, a  —  4(1 
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F 


X 


AT2 


va-  vc-  Eo,+T 


dTO; 

dT 


*o2 
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(A  14) 


bo2(t  -h)  =  -io2(t  -  h)  +  ^  go2,kk(t  -  h), 

k 

k  =  Da,  Dc,  T,  T.  (A  15) 
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Using  (A10)  and  (A13)  in  (A6),  it  can  be  found  the  equa¬ 
tions  for  the  anode  and  cathode  are. 


4(0 

4(0 


ga,a  ga,c  £a,T  £a,P 
gc,a  gc,c  gc,T  gc,P 

ba(t  -  h ) 
bc(t  -  h ) 


/«a(0\ 
Vc(t) 
Tit) 

V  Pit) ) 


(A  16) 


where, 

£a,a  — 
ga,c  — 
ga,  T  = 
&a,P  = 


gc,c  : 
gc,C  '■ 
~gc,T 
~gc,P 


gNi,  a  +  g02,a 
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:  gNi, T  +  g02,T 
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(A  17) 


baft  —  h)  =  —  Z?c(7  —  h)  =  b^it  —  h)  +  *o2(t  —  *).  (A  18) 
A.  3.  Thermal  terminal 

Construction  of  the  thermal  terminal  is  based  on  Eq.  (16), 
which  can  be  rearranged,  after  time-discretizing,  as: 


hit)  = 


cpm 

2*7X0 

r(0 

io2(0 

r(0 
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ar 

3£’o2 

ar 


(0 


(A  19) 


Thus,  we  obtain  the  RC  equation  for  the  thermal  terminal 


as: 


hit)  =  ( #T.a  <?T,c  <?T,T  <?T,p) 
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where, 
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(A24) 
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A.4.  Pressure  signal  terminal 

The  pressure  calculation  follows  Eqs.  (17)— (19).  A  time- 
discretized  form  of  the  pressure  is  given  as, 


pit)  = 


RTit ) 


3  V, 


8  L 


h 

—iNiit)  +  4^h2(^  -  h)  -  nn2it  -  2h) 


(A28) 


The  equation  for  the  pressure  signal  terminal  follows, 


Pit) 

1 p(0  =  — - pit). 

An 


(A29) 


where  p  denotes  the  pressure,  while  P  signifies  the  across 
variable.  Rp  is  the  output  impedance  of  the  pressure  terminal. 
To  accurately  monitor  the  pressure,  it  requires  /p(t)  ~  0  and 
Pit)  ~ pit).  Thus,  Rp  =  1.  The  RC  equation  for  this  terminal 
can  be  derived  as, 


zp(0=(gp,a  gp,c  gP,T  gp,p) 


^Pa(0\ 
Veit) 
Tit) 

V  PH)  1 


—  bvit  —  h) 


(A30) 
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where. 


gP,a  = 
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i>p(?  —  h)  =  —ip(t  —  h) 

+  f2gp,kk(t  -h),  k  =  va,  Vc,  T,  P. 


(A32) 


Combining  Eq.  (A  16),  (A20)  and  (A30)  yield  Eq.  (22). 
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